Anisotropy in Rupture Lines of Paper Sheets 
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We have experimentally investigated the fracture process in paper samples sub- 
mitted to a uniaxial force. Five types of paper sheets (newsprint, towel, sulfite, silk 
and couche papers) were fractured along two orthogonal orientations. In order to 
characterize the rupture lines of the paper sheets we utilized the Hurst exponent. 
Our results indicate a dependence of the Hurst exponent on the orientation of the 
paper sheets for samples of newsprint and, probably, towel and silk papers. For the 
other types of paper the Hurst exponent does not depend on the direction of crack 
propagation. 

PACS numbers: 62.20.Mk, 64.60.Fe, 05.40.-a 

I. INTRODUCTION of view of statistical mechanics, that utilize 

concepts such as percolation, fractals, scaling 
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rracture processes m disordered materials 
is a subject of intensive research and has at- Over the last decade fibrous materials 
tracted much scientific and industrial inter- have been attracting much attention. Several 



est 



ODD 



In any rupture experiment the models were elaborated to describe t 



le frac- 



presence of disorder is fundamental in the dy- ture process in these materials 

namics of crack propagation. It is naturally |ll|, Q Q Q Q Q- However, the small 

present in all materials and comes from a va- number of experimental data does not allow 

riety of different sources. Great theoretical a reliable comparison with the theoretical re- 

and experimental efforts have been done try- suits. 

ing to understand the process of crack forma- 

Experiments have shown that the crack 

tion and propagation in disordered materials. 

line in disordered materials can often be de- 



In particular, a significant fraction of these 

scribed by self-affine scaling |4 Il7l. Il8|. 

studies have been performed from the point 
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surface irregularity is often characterized by 
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rms value of the fluctuations of the surface 
height hi = h{xj) over a length scale e: 

i N—e -r i+e 

me) = -^E \ ^— r E (M*i) - <^» 2 

i=e \ ' j=i—e 

(1) 

where (h)i = Sjii- £ ^(%) is the average 
height around the position i. 

For a self-affine surface, the function h(x) 
is statistically invariant under an anisotropic 
scale transformation. This means that 
h(x) has the same statistical properties as 
k~ H h(kx), where H is known as the Hurst 
(roughness) exponent, which satisfies < 
H < 1. For self-affine structures, the rough- 
ness W over a range e satisfies the scaling 
law 

W(e) ~ e H . (2) 

The exponent H does not indicate how rough 
a surface is. It is a parameter that charac- 
terizes how the roughness, or the variance in 
the height, depends on the lateral scale over 
which it is measured. 

Some experimental works have claimed 
that the Hurst exponent H does not depend 
on the direction of crack propagation and as- 
sumes a universal value of 0.8 till 2o| - 
However, this universality was first ques- 
tioned by Milman et al. j^lj, who experi- 
mentally found with atomic force microscope 
measurements of fracture surface in crys- 
talline metals a Hurst exponent close to 0.5, 
a value significantly smaller than 0.8. An ex- 



ponent H = 0.5 was also observed in mate- 
rials displaying intergranular fractures, such 
as sandstone [22]. From the theoretical point 
of view, numerical models have been elabo- 
rated in order to evaluate the Hurst exponent 
H . Computer simulations have shown that 



2J and that 



H ~ 0.7 in two dimensions 

0,1 

H r ang es from 0.4 to 0.5 in three dimensions 
I25I \l\ . Nowadays there is a conjecture relat- 
ing the value of the Hurst exponent H to the 
speed of crack propagation through the sam- 
ple |26j]. The higher value (H = 0.8) has been 
associated with a high speed of crack propa- 
gation and interpreted dynamic regime. 
In contrast, the smaller value (H = 0.5) was 
related with a quasi-static regime, in which 
the dynamic effects of crack propagation are 
negligible. Experimentally, the smaller and 
the higher values of H are connected to the 
eng th scale at which the crack is examined 
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Parisi and Caldarelli [29( have studied the 
fracture surface of three dimensional sam- 
ples using a model for quasi-static fractures 
known as Born Model. They have found a 
Hurst exponent H = 0.5, which does not de- 
pend on the direction of crack propagation. 

An experimental procedure used to inves- 
tigate the fracture process of fibrous materi- 
als is the rupture of a paper sheet into two 



parts 
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Paper is a good example 



of a disordered material due to its structural 
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non-uniformity, which can be confirmed even 
with the naked eye. When light is transmit- 
ted through it, patches of lighter and darker 
regions can be seen. Such darker and lighter 
parts correspond, respectively, to higher and 
lower local densities. Paper is made from 
cellulose fibers found in wood, the most fre- 
quently used being eucalyptus and pine. Dif- 
ferent types of wood can be utilized to pro- 
duce cellulose fibers with specific character- 
istics, therefore producing papers with dif- 
ferent proprieties. Generally, the fiber net- 
work structure is highly disordered since their 
fibers are randomly positioned and oriented. 
Nevertheless, in some types of paper, the fab- 
rication process is such that the fibers acquire 
some order, having a tendency to align along 
one direction. The fibers are held together 
by their ability to form hydrogen bonds with 
each other. The paper strength is a funda- 
mental feature belonging to that group of 
physical properties associated with the pa- 
per manufacturing process. Indeed, paper 
strength is influenced by the kind, quality, 
and treatment of the constituent fibers and 
by the way the sheet has been formed on the 
paper machine. 

The analysis of rupture in paper samples 
is very suitable because the best method to 
calculate the Hurst exponent demands the 
knowledge of the profile, i.e, of the func- 
tion h = f(x). In the case of a two- 



dimensional fracture such a height profile is 
easily obtained. However, the evaluation of 
the Hurst exponent shows a significant dis- 
persion, therefore demanding the analysis of 
a great number of samples. There are in the 
literature some experiments that investigated 
the rupture process in paper samples. 

Maunuksela et al. [32], performed exper- 
iments on the dynamics and kinetic rough- 
ening of one-dimensional slow-combustion 
fronts using three types of paper. They con- 
cluded that well-controlled experiments lead 
to a clear asymptotic scaling that unequiv- 
ocally lies in the thermal KPZ universality 
class. 

Kertsz et al. Q| have done experiments 
using a tensile testing machine to study the 
morphology of rupture lines in paper. The 
samples were submitted to traction and the 
velocity of crack propagation was kept con- 
stant until the breaking process separated the 
sample into two pieces. They have studied 
five different kinds of paper and found that 
the Hurst exponent characterizing the rup- 
ture lines did not depend significantly on the 
type of paper used in the rupture process. 
The results obtained indicated a Hurst expo- 
nent ranging from 0.63 to 0.72, which did not 
depend on the direction of crack propagation. 

Salminen et al. [33] have studied frac- 
ture in paper via acoustic emission analysis. 
In this case, the release of acoustic energy 
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is related to irreversible deformation, micro- II. EXPERIMENTAL PROCEDURE 



craks, and, perhaps, to plasticity. It has 

been verified that the acoustic emission of We have carried out rupture experi- 

energy and the waiting times between acous- ments using five different types of commer- 

tic events followed power-law distributions, cially available papers: newsprint ((36.8 ± 

This remained true while the strain rate was 0.2) g/m 2 ), towel ((31.7 ± 0.1) g/m 2 ), 

varied by more than 2 orders of mag nitude. silk ((51.8 ± 0.5) g/m 2 ), couche ((82.7 ± 

Recently, Salminen et al. 31] have analyzed 0.5) g/m 2 ) and sulfite ((78.1 ± 0.4) g/m 2 ). 

6.5 m long fracture lines of paper as an exam- Two directions, x and y, perpendicular to 

pie of crack propagation involving disorder, each other, were assigned to the paper, along 

They have observed that the roughness ex- the width and the length of the sheets respec- 

hibits a power-law scaling, with an exponent tively. We refer to the x as the horizontal 

close to 0.60 ±0.1 for all samples. direction and to the y as the vertical direc- 
tion (Fig. 1). Square samples of 20 x 20 cm 
were cut from the paper sheets and the orien- 

In this paper we have studied the depen- tation of the directions x and y were marked 

dence of the rupture line morphologies on the on them - Each sample was fixed at both ends 

direction of crack propagation. Several sam- b y two parallel clamps, either with their x or 

pies of five different types of paper were sub- V direction parallel to the clamps. The up- 

mitted to a constant load until their com- P er clam P was maintained fixed and in the 

plete rupture. Two directions, perpendicular lower one a static load F was applied. The 

to each other, were chosen for crack propaga- same load F ~ 200 N was used in all exper- 

tion and the fracture surfaces obtained were iments. In order to initiate the propagation 

characterized by their Hurst exponent. of the crack a notch 2 cm lon S was cut into 

the paper in the middle of one of the sides. 

After fractured, the samples were digitized 

The paper is organized as follows. The in a scanner of 450 dpi resolution with 256 

next section is devoted to present the exper- gray levels and a black background. During 

imental procedure used in the rupture pro- the scanning process the beginning and the 

cess. In section III, the obtained results are end of the rupture line were discarded. The 

presented and analyzed and, finally, the last gray level distribution near the front of the 

section is devoted to conclusions. digitized figure ranges from about 170 on the 
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FIG. 1: Schematic representation of the experi- 
mental apparatus used in the rupture of paper. 
The direction of crack propagation is indicated 
in each figure, a) horizontal direction and b) 
vertical direction. 

side belonging to the paper to about 70 out- 
side of the paper (on a scale that ranges from 
to 256). The number of colors used to dig- 
itize the samples was later decreased to only 
2 (black and white) using the function black 
and white from CorelPhoto-Paint 11. The 
threshold was adjusted until just the smaller 
fibers were deleted from the front. The digi- 
tized patterns were analyzed by a home made 
computer program which allowed us to obtain 
the functions h(x) which represent the pat- 
terns of the investigated rupture lines. Fig- 
ure 2 shows typical functions h(x) obtained 
for a newsprint sample along the horizontal 
and vertical directions. 

III. EXPERIMENTAL RESULTS 

From the digitized patterns we evaluated 
the roughness W (using Eq. (1)) of the rup- 
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FIG. 2: Functions h(x) for two samples of 
newsprint: (a) horizontal and (b) vertical direc- 
tion. The horizontal and vertical scales are in 
pixels. 

ture lines. A log-log plot of W versus e 
presents a region where the data have a linear 
behavior over at least two decades. The slope 
of this region is equal to the Hurst exponent. 

Figure 3 shows typical results for samples 
of newsprint, silk and towel papers along the 
two perpendicular directions. The slopes of 
the two straight lines shown are equal to 0.60 
and 0.72 for the newsprint, 0.85 and 0.94 for 
the silk paper and 0.71 and 0.78 for the towel 
paper. The lower value corresponds to the 
value of the Hurst exponent along the hor- 
izontal direction and the higher one to the 
value of the Hurst exponent along the verti- 
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FIG. 3: Roughness of rupture lines for two sam- 
ples of: (a) newsprint, (b) silk and (c) towel pa- 
pers. The rupture lines are oriented along two 
different directions: vertical (circles) and hori- 
zontal (up triangles). The slopes of the linear 
parts allow us to evaluate the Hurst exponents. 

Table I presents the values obtained for 
the Hurst exponents calculated for each type 
of paper, along the horizontal and vertical di- 
rections. For all five types of paper, the data 
were averaged over 30 statistically indepen- 
dent samples. 

Our results indicate a dependence of the 
Hurst exponent upon the direction of crack 



Paper H y H x 

newsprint 0.73 ± 0.02 0.65 ± 0.02 
couche paper 0.74 ± 0.04 0.74 ± 0.03 
towel paper 0.78 ± 0.02 0.71 ± 0.03 
sulfite paper 0.80 ± 0.04 0.77 ± 0.04 
silk paper 0.94 ± 0.02 0.90 ± 0.01 

TABLE I: Hurst exponent for different kinds of 
paper sheet obtained along two directions (H y = 
Hurst exponent for the vertical direction and 
H x = Hurst exponent for the horizontal direc- 
tion). 

propagation for the samples of newsprint, 
towel and, probably, silk papers. This de- 
pendence is related to the inhomogeneity of 
the paper. For these three types of paper the 
Hurst exponent is larger along the vertical di- 
rection than along the horizontal one. This 
indicates that in these papers, the propaga- 
tion speed of a crack should be larger along 
the vertical direction. 

For the other two types of paper it was 
observed that, within the experimental er- 
ror, the Hurst exponent of the studied sam- 
ples did not depend on the direction of crack 
propagation. Notice that in Table I the val- 
ues obtained for the Hurst exponents along 
both directions are higher than 0.5. This re- 
sult indicates that the fracture process inves- 
tigated here can be identified as belonging to 
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the dynamic regime and, also, corresponds to In the method of the best linear least- 



measurements at large length scales 



square fitting, the roughness W(e) over a 



The results in Table I, for the newsprint, lenght scale e is given by 
clearly reveal that the Hurst exponent de- ^ n 

pends on the direction along which it is mea- ^( e ) jyX^^ 6 ) CO 



and the local roughness u>i(e) is defined as 



sured. In order to confirm this result we 
have increased the scanner resolution to 800 

dpi and used three other test methods for u> 4 2 (e) = jz rr X! \hj ~ ( a i( e ) x j + ^( e ))f- 

the crack roughness analysis: max-min dif- (8) 

ference |34j, power spectrum [35| and best a .( e ) an d fc.( e ) are the linear fitting coeffi- 

linear least-square fitting |36j. c i ents to the displacement ratio data on the 

The max-min method consists of the com- interval [i-e,i + e] centered at the position i. 

putation of h max (e) } defined as the difference For a self-affine profile, the roughness W of 

between the lowest and the highest heights h the fracture surface satisfies the scaling law 



W(e) ~ e H . (9) 



within a certain window of size e, averaged 
over all possible origins x of the window. For 
a self-affine profile, h max scales as a function The valueg of the Hurgt exponentS) es _ 

of e according to timated using the four methods, for the 

^ ^ e H ^ newsprint are given in Table II. Notice that 

the three methods agree on the values of the 

The power spectrum method measures the TT , , ^ . £ , . , 

Hurst exponents. 1 his fact indicates that the 

power spectrum of the interface, and not the , r , ir ^ i 

crack surfaces are truly selt-attine and con- 
height h(x). In this case it is calculated the n ,i • , , , , , , TT 

v ' firms the conjecture that the Hurst exponent 

structure factor, or power spectrum , , rr 

depends on the direction of fracture. 

(jfc) = (h(k)h{-k)) , (4) Figures 4(a) and 4(b) show optical mi- 

crophotographs of newsprint and couche pa- 
per samples, respectively. For the newsprint 
h(k) = — y^\h(x) — (h)i)]exp(ik.x) , (5) it can be clearly seen that there are many 



where 

= L 1 / 2 



X 



long fibers running, preferentially, along the 
horizontal direction of the figure. This ori- 
entation coincides with the horizontal direc- 
s(k) ~ k~( 2H+1 \ (6) tion of the paper. For the couche paper the 



and k is the wave factor. The structure factor 
scales, for a self-affine profile, as 
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process can occur due to two different pos- 

method H y H x 

sitrte fracture modes: tensile fracture of the 
rms 0.69 ± 0.02 0.61 ± 0.02 

hydrogen bonds between the fibers and shear 
max-min 0.73 ± 0.02 0.63 ± 0.02 

fracture of the fibers. When the newsprint 
power spectrum 0.71 ± 0.03 0.63 ± 0.03 

is fractured along the vertical direction the 

best linear least-square fitting 0.68 ± 0.03 0.58 ± 0.03 

crack has to propagate in a direction perpen- 
dicular to the fibers. In this case, the hy- 

TABLE II: Hurst exponent obtained for 

drogen bonds between the fibers are mainly 

newsprint using different methods. 

subjected to shear fracture mode while the 
fibers are also long but they do not appear fibers are mainl y subjected to tensile frac- 
to be oriented along any preferential direc- ture mode - A S ain the ru P ture of ^dividual 
tion, being more intertwined than those of fibers is not necessary. If the hydrogen bonds 
the newsprint. The orientation of the fibers between them break, they can simply slide 
is produced during the paper fabrication pro- P ast each other - The separation of the fibers 
cess and seems to play an important role in 1S P robab ly more difficult m this process as 
the roughness of the rupture lines. The ac- compared to that for the horizontal direction 
tual mechanism of fracture is not yet com- since > as the fibers slide P ast each other > h y~ 
pletely understood at the microscopic level. dr °g en bonds that were broken can be formed 
For the paper to be fractured, the rupture of a § ain - Rupture of individual fibers should be 
the hydrogen bonds between the fibers must even more difficult since the force necessary 
occur in order that they can be separated to rupture one fiber must be greater than that 
from one another. However, the rupture to break a h J^gen bond. Consequently, the 
of individual fibers seems to be not neces- crack has to chan S e its direction of propaga- 
sary. When the newsprint is fractured along tion man y times > tr y in S to find a weaker P oint 
the horizontal direction the crack can prop- to fracture and the m P ture line formed is 
agate mainly through the separation of one rou S her than that formed alon S the horizon- 
fiber from another by breaking the hydrogen tal direction. Microscopic examination of the 
bonds. The rupture of individual fibers, if it ru P ture lines show that there are much more 
occurs at all, does not have to happen very lon S fibers protruding from the paper when 
often since most of the fibers are aligned par- the P a P er is fractured along the vertical direc- 
allel to the direction of crack propagation. In tion than when [t is fractured along the hori- 
this case, we can suppose that the rupture zontal direction. This indicates that, indeed, 
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the fibers are just being pulled apart without 
being fractured. Unfortunately, comparison 
of the end of the protruding fibers with the 
ends of the fibers away from the rupture line 
could not reveal if the former where broken 
during the fracture process. It is possible to 
make an analogy between paper fracture and 
microscopic failure mechanisms of a unidirec- 
tional composite material Q if the hydrogen 
bonds are considered as being the matrix. 



For the other types of paper that pre- 
sented a dependence of the Hurst exponent 
on the orientation of the paper sheet, it was 
also possible to see the same kind of preferen- 
tial fiber orientation, although much less pro- 
nounced than for the newsprint. This, proba- 
bly, is the reason why the difference between 
H x and H y for the towel and the silk papers 
is smaller than for the newsprint. For the two 
types of paper that did not show dependence 
between the Hurst exponent and the orien- 
tation of the paper sheet it was not possible 
to detect any orientation of the fibers. There 
seems to be, also, a relation between the size 
of the fibers and the value of the Hurst ex- 
ponent. Papers composed of long fibers have 
smaller exponents than those of papers com- 
posed of short fibers. 



FIG. 4: Optical microphotographs for: (a) 
newsprint and (b) couche paper. Both papers 
are oriented with their y direction along the 
length of the page. 

IV. CONCLUSION 

In this paper we have experimentally in- 
vestigated the fracture process in five dif- 
ferent types of paper submitted to a static 
load. The rupture line of the samples were 
characterized by the Hurst exponent, calcu- 
lated along two orthogonal orientations of the 
paper: vertical and horizontal. This proce- 
dure was utilized to investigate the effect of 
anisotropy on the Hurst exponent. 

For the samples of newsprint our results 
clearly indicate that the Hurst exponent de- 
pends on the orientation of the paper. This 
dependence is associated with the orientation 
of the fibers within the paper, which is de- 
termined by its fabrication process. For the 
samples of silk and towel papers the differ- 
ence between the Hurst exponents is small, 
almost within the experimental error, indi- 
cating a probable dependence of the Hurst 
exponent upon the orientation of the paper. 

For all paper samples studied the rupture 
process is very rapid. This can be verified 
by the value of the Hurst exponent, higher 
than 0.5 for all samples tested. Therefore, 
the fracture process can be classified as a dy- 
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namic one. 

Contrary to the results reported on the lit- 
erature, our results support the idea that the 
value of the Hurst exponent for paper sheets 
is not universal. It depends on both the type 
of paper used and, more interestingly, it can 
also depend on the direction along which it 
is measured for the same paper sheet. 
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